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Abstract 

A measurement of the inclusive Z — )■ rr cross-section in pp collisions at ^/s = 7 TeV 
is presented based on a dataset of 1.0 fb~^ collected by the LHCb detector. Can- 
didates for Z — 7- rr decays are identified through reconstructed final states with 
two muons, a muon and an electron, a muon and a hadron, or an electron and 
a hadron. The production cross-section for Z bosons, with invariant mass be- 
tween 60 and 120 GeV/c^, which decay to r leptons with transverse momenta 
greater than 20 GeV/c and pseudorapidities between 2.0 and 4.5, is measured to 
be (Jpp^z^TT = 71.4 lb 3.5 ib 2.8 it 2.5 pb; the first uncertainty is statistical, the 
second is systematic, and the third is due to the uncertainty on the integrated lu- 
minosity. The ratio of the cross-sections for Z —)• rr to Z —t- ^u/i is determined to 
be 0.93 ± 0.09, where the uncertainty is the combination of statistical, systematic, 
and luminosity uncertainties of the two measurements. 
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1 Introduction 



The measurement of the production cross-section for Z bosons* in proton-proton i^p'p) col- 
hsions constitutes an important verification of Standard Model predictions. Since lepton 
universahty in Z decays has been tested to better than 1% at LEP [1], any deviation ob- 
served at the LHC would be evidence for additional physics effects producing final state 
leptons. In particular, r-lepton pairs can be important signatures for supersymmetry, 
extra gauge bosons, or extra dimensions [2, 3, 4]. The LHCb experiment has previously 
measured the cross-section for Z — )• /x/x [5] with both leptons having transverse momen- 
tum (pt) above 20 GeV/c and an invariant mass between 60 and 120 GcV/c^. Here a 
complementary measurement in the decay mode Z —)■ rr is presented. This measurement 
extends the Z — >■ rr cross-section measurements from the central pseudorapidity range 
covered by ATLAS (|?7| < 2.4) [6] and CMS (I77I < 2.3) [7] into the forward region covered 
by the LHCb experiment (2 < 77 < 4.5). 

2 Detector and datasets 

The LHCb detector [8] is a single-arm forward spectrometer designed for the study of 
particles containing 6 or c quarks. The detector includes a high precision tracking system 
consisting of a silicon-strip vertex detector (VELO) surrounding the pp interaction re- 
gion, a large-area silicon-strip detector (TT) located upstream of a dipole magnet with a 
bending power of about 4 Tm, and three stations of silicon-strip detectors (IT) and straw 
drift tubes (OT) placed downstream. The combined tracking system has a momentum 
resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and an impact 
parameter resolution of 20 //m for tracks with high px. 

Charged hadrons are identified using two ring-imaging Cherenkov detectors. Pho- 
ton, electron and hadron candidates arc identified by a calorimeter system consisting of 
scintillating-pad (SPD) and pre-shower detectors (PRS), an electromagnetic calorimeter 
(ECAL) and a hadronic calorimeter (HCAL). Muons are identified by a system composed 
of alternating layers of iron and multiwire proportional chambers. The trigger consists of 
a hardware stage, based on information from the calorimeter and muon systems, followed 
by a software stage that applies a full event reconstruction. The hardware stage imposes 
a global event requirement (GEC) on the hit multiplicities of most sub-detectors used in 
the pattern recognition algorithms to avoid overloading of the software trigger by high 
occupancy events. 

This analysis uses data, corresponding to an integrated luminosity of 1028 ± 36 pb~^, 
taken at a centre-of-mass energy of 7 TeV. The absolute luminosity scale was measured 
periodically throughout the data taking period using Van der Meer scans [9] where the 
beam profile is determined by moving the beams transversely across one another. A 
beam-gas imaging method was also used where the beam profile is determined through 
reconstructing beam-gas interaction vertices near the beam crossing point [10]. Both 

*Here, the Z is used to indicate production from Z bosons, photons, and their interference. 
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methods provide similar results and the integrated luminosity is determined from the 
average of the two, with an estimated systematic uncertainty of 3.5% [11]. The primary 
systematic uncertainty of 2.7% is due to the beam current measurement, shared between 
the two methods. 

Simulated data samples are used to develop the event selection, determine efficien- 
cies, and estimate systematic uncertainties. Each sample was generated using an LHCb 
configuration [12] of Pythia 6.4 [13] with the CTEQ6L1 leading-order PDF set [14] and 
passed through a Geant4 [15] based simulation of the LHCb detector [16]. Trigger 
emulation and full event reconstruction were performed using the LHCb reconstruction 
software [17]. Additional samples, without detector simulation or event reconstruction, 
are used to study the signal acceptance and were generated using Pythia 8.1.55 [18], 
Herwig-F-F 2.5.1 [19], and Herwig-F-F with the POWHEG method [20]. 

3 Event selection 

The signatures for Z — )■ rr decays considered in this analysis are two oppositely-charged 
tracks, consistent with an electron, muon, or hadron hypothesis, having large impact 
parameters with respect to the primary vertex of the event. 

Tracks are reconstructed in the VELO and extrapolated to the IT/OT sub-detectors; 
any TT sub-detector hits consistent with the track are added and a full track fit is per- 
formed. Only tracks with fit probabilities greater than 0.001 are considered. 

Tracks are extrapolated to the calorimeters and matched with calorimeter clusters. 
Electron candidates arc required to have a PRS energy greater than 0.05 GeV, a ratio 
of ECAL energy to candidate momentum, E/pc, greater than 0.1, and a ratio of HCAL 
energy to candidate momentum less than 0.05. Any electron candidate momentum is 
corrected using bremsstrahlung photon recovery [21]. Since the ECAL is designed to 
register particles from 6-hadron decays, calorimeter cells with transverse energy above 
10 GcV saturate the electronics, and lead to degradation in the electron energy resolution. 

Hadron candidates are identified by requiring the ratio of HCAL energy to track 
momentum to be greater than 0.05. Due to the hmited HCAL acceptance, the candidate 
track is required to have a pseudorapidity of 2.25 < r] < 3.75. 

Muon candidates are identified by extrapolating tracks to the muon system down- 
stream of the calorimeters and matching them with compatible hits. Muon candidates 
are required to have a hit in each of the four stations and consequently will have traversed 
over 20 hadronic interaction lengths of material. 

The data have been collected using two triggers: a trigger which selects muon candi- 
dates with a pt greater than 10 GeV/c; and a trigger which selects electron candidates 
with pt greater than 15 GeV/ c. 

The analysis is divided into five streams, labelled t^t^, TuTc TsT^, T/^Th, and TgTh, 
defined such that the streams are exclusive. The first r lepton decay product candidate 
is required to have pt > 20 GeV/ c and the second is required to have px > 5 GeV/ c. The 
following additional kinematic and particle identification requirements are specific to each 
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analysis stream: 

• Tfj,Tn requires two oppositely-charged muons where at least one triggered the event. 
The muon with the larger is considered as the first r lepton decay product 
candidate. 

• r^Te requires a muon that triggered the event and an oppositely-charged electron. 

• TeT^ requires an electron and an oppositely-charged muon with < 20 GeV/c. 
Either lepton can trigger the event. 

• TixTh requires a muon that triggered the event and an oppositely-charged hadron. 

• TgTh requires an electron that triggered the event and an oppositely-charged hadron. 

In pp collisions the cross-section for hadronic QCD processes is very large. These 
events can pass the above requirements either due to semileptonic c- or 6-hadron decays 
or through the misidentification of hadrons as leptons. 

Signal decays, coming from an on-shell Z, tend to have back-to-back isolated tracks 
in the transverse plane with a higher invariant mass than tracks in QCD events. The 
absolute difference in azimuthal angle of the two r lepton decay product candidates, 
|A$|, is required to be greater than 2.7 radians and their invariant mass is required to be 
greater than 20 GeV/c^. 

Tracks in QCD events also tend to be associated with jet activity, in contrast to signal 
events where they are isolated. An isolation variable, Jp^, is defined as the transverse 
component of the vectorial sum of all track momenta that satisfy \/A0^~+~A?p < 0.5, 
where A0 and Ar) are the differences in and r) between the r lepton decay product 
candidate and the track. The track of the r lepton decay product candidate is excluded 
from the sum. Both r lepton decay product candidates are required to have /p^ < 2 GeV/ c 
for the T^T^, T^Te, and Tgr^ analysis streams and Ip^ < 1 GeV/c for r^r/j and TeTh due to 
the larger QCD backgrounds. 

The hfetime of the r lepton is used to separate signal from prompt backgrounds. The 
signed impact parameter for a track is defined as the magnitude of the track vector of 
closest approach to the primary vertex signed by the ^-component of the cross product 
between this vector and the track momentum. The impact parameter significance, IPS, 
is then defined as the absolute sum of the signed impact parameters of the two r lepton 
decay product candidates, divided by their combined uncertainty. The IPS is required to 
be greater than 9 for the t^t^, T^Th, and T^Th analysis streams while no IPS requirement 
is placed on the r^Tg or Tgr^ streams. 

In the r^r^ analysis stream an additional background component arises from Z — )■ nn 
events. This produces two muons with similar px, most of which also have an invariant 
mass close to the Z mass. In contrast, signal events tend to have unbalanced px and a lower 
invariant mass due to unreconstructed energy from neutrinos and neutral hadrons. The px 
asymmetry, Ap^, is defined as the absolute difference between the px of the two candidates 
divided by their sum. For the r^r^ analysis stream the Ap^ is required to be greater than 
0.3 and the di-muon invariant mass must lie outside the range 80 < M^^ < 100 GeV/c^. 
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4 Background estimation 



The invariant mass distributions for the selected Z ^ tt candidates, the simulated signal, 
and the estimated backgrounds for the five analysis streams are shown in Fig. 1, where no 
candidates are observed with a mass above 120 GeV/c^. Five types of background have 
been considered: generic QCD; electroweak, where a high pt lepton is produced by a 
or Z boson and the second candidate r lepton decay product is misidentified from the 
underlying event; tt, where two hard leptons are produced from top decays; WW , where 
each W decays to a lepton; and Z II for the t^t^, t^t/j, and r^Th streams, where for the 
T^Th stream a single muon is misidentified as a hadron and for the T(,Tf, stream an electron 
is misidentified. The ti and WW backgrounds are estimated using simulation and found 
to be small for all final states. 

The QCD and electroweak backgrounds are estimated from data. A signal-depleted 
control sample is created by applying all selection criteria but requiring that the r lepton 
decay product candidates have the same-sign (SS) charge. The QCD and electroweak 
background events in this sample, N^^) ^i^d iV|wK) obtained by fitting template 
shapes to the distribution of the difference between the of the first and second r 
lepton decay product candidates. The template shape for the electroweak background is 
taken from simulation. To determine the shape of the QCD contribution, the isolation 
requirement is reversed such that Jp^ > 10 GeV/c. The number of candidates for each 
background category in the signal sample is then calculated as A^qcd = /qcdA/'qcd ^^^1 
Newk — /ewk-/Vewk) where /qcd and /ewk are the ratio of opposite-sign to same-sign 
events for QCD and electroweak events respectively. Both /qcd and /ewk are determined 
as the ratio of opposite-sign to same-sign events satisfying the template requirements. The 
uncertainties on the QCD and electroweak backgrounds are estimated by combining the 
statistical uncertainty on the fraction with the uncertainties from the fit used to determine 
N^CB and iV||,K- 

The number of Z — >■ /x/i background events for the r^r^ stream is obtained by applying 
all selection criteria except for the 80 < M^^ < 100 GcV/c^ requirement. This produces a 
sample with a clear peak around the Z mass as shown in Fig. 1(a). A template for Z — )■ fifi 
events is obtained from data by applying the event selection, but requiring prompt events 
with IPS < 1. The template is normalised to the number of events within the t^uT^ sample 
with IPS > 9 and within the invariant mass range 80 < M^^ < 100 GeV/c^. The Z ^ jiji 
background is the number of events in the normalised template outside this mass range. 
The uncertainty on this background is estimated from the statistical uncertainty on the 
normalisation factor. 

The Z ^ nil process also contributes a small background to the T^Th stream when 
one of the muons is misidentified as a hadron. This is evaluated by applying the t^t^ 
selection but requiring a second identified muon rather than a hadron, and scaling this 
by the probability for a muon to be misidentified as a hadron. The latter is found from a 
sample of Z ^ ii/i events that have been selected by requiring a single well defined muon 
and a second isolated track, which give an invariant mass between 80 and 100 GeV/c^; 
(0.06 ± 0.01)% of these tracks pass the hadron identification requirement. 
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Figure 1: Invariant mass distributions for the (a) t^t^, (b) r^Tg, (c) Tgr^, (d) r^Th, and (e) 
TeTh candidates with the excluded mass range indicated for r^r^. The Z — j- rr simulation 
(solid red) is normalised to the number of signal events. The QCD (horizontal green), 
electroweak (vertical blue), and Z (solid cyan) backgrounds are estimated from data. 
The ti (vertical orange) and WW (horizontal magenta) backgrounds are estimated from 
simulation and generally not visible. 
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Table 1: Acceptance factors, branching fractions, selection efficiencies, numbers of back- 
ground and observed events for each Z ^ tt analysis stream. 



Stream 


A 


B [%] 




A^bkg 


N 




0.405 ± 0.006 


3.031 ±0.014 


0.138 ±0.006 


41.6 ±8.5 


124 




0.248 ± 0.004 


6.208 ±0.020 


0.517 ±0.012 


129.7 ±4.9 


421 


TeTfj. 


0.152 ±0.002 


6.208 ±0.020 


0.344 ±0.016 


56.6 ±3.3 


155 




0.182 ±0.002 


16.933 ±0.056 


0.135 ±0.004 


53.3 ±0.8 


189 


TeTh 


0.180 ±0.002 


17.341 ± 0.057 


0.082 ± 0.004 


36.6 ±0.9 


101 



Similarly, a small Z — )• ee background can contribute to the rgT/j stream when one of 
the electrons is misidentified as a hadron. This is evaluated by applying the r^Th selection 
but requiring a second identified electron rather than a hadron, and scaling this by the 
probability for an electron to be misidentified as a hadron. The electron mis-identification 
is found from simulated Z — > ee events to be (0.63 ± 0.02)%. 

5 Cross-section measurement 

The pp Z ^ TT cross-section is calculated within the kinematic region 60 < M^^ < 
120 GeV/c2, 2.0 < r]^ < 4.5, and > 20 GeV/c using 

where N is the number of observed candidates and iVj^^g is the estimated background from 
source j. The integrated luminosity is given by A is an acceptance and final state 
radiation correction factor, B is the branching fraction for the r-lepton pair to decay to 
the final state, and £sei is the selection efficiency. A summary of these values for each final 
state is given in Table 1. The reconstruction efficiency, £rec; is calculated using simulation 
or data for each event, assuming that it is signal, and depends on the momentum and 
pseudorapidity of the r lepton decay product candidates. (l/£rec)"' indicates the average 
value of l/^rcc for background source j. 

The integrated luminosity of the datasets for the t^t^, T^Tg, and r^r/j samples is 1028 ± 
36 pb""*^, while the TgT^ and TeTh final state datasets have an integrated luminosity of 
955 ± 33 ph'\ 

5.1 Acceptances and branching fractions 

The acceptance factor. A, is used to correct the kinematics of each analysis stream to 
the kinematic region 60 < M^r < 120 GeV/c^ 2.0 < 77^ < 4.5, and > 20 GeV/c. This 
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region corresponds to the detector fiducial acceptance and allows a comparison with the 
LHCb Z fj^fj, measurement [5]. The acceptance factor is taken from simulation and is 
defined as the number of Z —>■ rr events where the generated r lepton decay products 
fulfil the kinematic requirements described in Sect. 3, divided by the number oi Z ^ tt 
events where the generated r leptons lie within the kinematic region defined above. 

For each final state the acceptance factors are calculated at leading-order using fully 
modelled hadronic decay currents and spin correlated r lepton decays with final state 
radiation in Pythia 8 and HERWIG++, and at next-to-leading-order using the POWHEG 
method implemented in HerwiG-F-|-. For Pythia 8 the CTEQ5L leading-order PDF 
set [22] was used, while for Herwig++ the MSTW08 PDF set [23] was used. The mean 
of the maximum and minimum values from the three generators is taken as the acceptance 
factor and is given in Table 1. The uncertainty is taken as half the difference between the 
maximum and minimum values. 

The branching fractions are calculated using the world averages [24] and are given 
in Table 1. The r lepton to single charged- hadr on branching fraction is the sum of all 
r lepton decays containing a single charged hadron. The final states presented in this 
analysis account for 44% of all expected Z — > rr decays. 

5.2 Selection efficiency 

The event selection efficiency, £sei, is the product of the efficiencies described below. Each 
efficiency is determined from either data, or simulation which has been calibrated using 
data. The resulting Sgei for each stream is given in Table 1. 

The kinematic efficiency, £km, is obtained from simulation and is the number of events 
fulfilling the kinematic requirements of Sect. 3 at both the simulated and reconstructed 
level divided by the number of events passing the requirements at the simulated level. 
The efficiency is consistent with unity for the r^r^ and t^t/j analysis streams. For streams 
involving electrons, £kin is significantly lower due to the saturation of the ECAL. This 
results in electrons being reconstructed with lower momenta than their true momenta 
due to incomplete bremsstrahlung recovery. In the r^Tg, TgT^, and r^Th streams, 5kin is 
(99.3 ± 1.0)%, (66.8 ± 1.9)%, and (67.0 ± 1.3)% respectively. The uncertainties come from 
the statistical uncertainty of the Z — > rr simulation and the cahbration of the electron 
momentum scale which has been obtained by comparing the px spectrum of Z — > ee 
events in data and simulation [25]. 

The efficiency of the isolation requirement, £/p^, for each analysis stream is taken 
from Z ^ TT simulation, and calibrated to data by multiplying si^^ by the ratio of the 
efficiency obtained in Z — )• /x/i data to Z — >■ r^r^ simulation. The systematic uncertainty 
on £/p^ is estimated as the difference between the efficiencies obtained from Z ^ ii/j, 
simulation and Z r^r^ simulation. 

The efficiency of the impact parameter significance requirement, eips, is evaluated 
from Z ^ TT simulation. A comparison of the IPS distributions in Z ^ ^fi events 
from data and simulation show that the impact parameter resolution is underestimated 
by (12 ± 1)% in simulation, and so the simulated Z — > rr events are corrected by this 
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factor. The systematic uncertainty on eips is determined by re-calculating the efficiency 
in Z TT simulation with the scale factor varied by its uncertainty. 

The efficiency of the azimuthal angle separation requirement, £|a*|, and px asymmetry 
efficiency requirement, EAp^, are evaluated from simulation. The systematic uncertainty 
on each is taken as the difference in the evaluation of these efficiencies in Z — )■ /i/x data 
and simulation, combined in quadrature with the statistical uncertainty from the Z ^ tt 
simulation. 

5.3 Reconstruction efficiency 

The reconstruction efficiency, £rec; is the product of the GEC, trigger, and tracking and 
identification efficiencies for both r lepton decay product candidates. The tracking effi- 
ciency is the probability for reconstructing the track and the identification efficiency is 
the probability for the track to be identified by the relevant sub-detectors. All efficiencies 
determined from data have been checked against simulation and found to agree within 
the percent level. 

The GEC efficiency, EgeC: is a correction for the loss due to the rejection by the 
hardware trigger of events with an SPD multiplicity of greater than 600 hits. For muon 
triggered events, the efficiency has been evaluated to be (95.5 ± 0.1)% from Z ^ jiji data 
events using a hardware di-muon trigger with a relaxed SPD requirement of 900 hits. For 
electron trigger events, the efficiency is estimated to be (95.1 ± 0.1)% by comparing the 
hit multiplicities in Z ^ /i/i and Z — >■ ee events. 

The muon and electron trigger efficiencies, £trg, are evaluated in bins of momentum 
using a tag-and-probc method on Z ^ ££ data events, which have been selected requiring 
two reconstructed and identified muon or electron candidates with an invariant mass 
consistent with that of the Z. In the events the triggered lepton is taken as the tag lepton, 
and the other as the probe lepton. The event topologies for Z ^ ii and Z — > rr events 
are nearly identical except for the momenta of the final state particles and so the trigger 
efficiency is calculated only as a function of the probe momentum below 500 GcV/c. The 
trigger efficiency is the fraction of events where the probe has also triggered, and varies as 
a function of probe momentum between 75% and 80% for the muon trigger and between 
62% and 75% for the electron trigger. The trigger efficiency uncertainty for each bin in 
momenta is taken as the statistical uncertainty. 

The tracking efficiency, Sti-k, is also evaluated for muons using a tag-and-probc method 
on the Z ^ /jfi data. The tag must satisfy all the muon reconstruction and identification 
requirements. The probe is reconstructed from a track segment in the muon chambers 
that has been associated to a hit in the TT sub-detector, which is not required in the track 
reconstruction. Events with a tag and probe mass consistent with the on-shell Z mass are 
used. The tracking efficiency is evaluated as the number of events with a reconstructed 
probe track over the total number of events. For lower px tracks, masses consistent with 
the J/ijj are used. 

The J/ip — )■ yU/x topology differs from the Z — >■ rr topology in both pseudorapidity 
and momentum, and so the J/ip muon tracking efficiencies are evaluated in bins of both 
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variables. The muon tracking efficiency is found to vary between 85% and 93%. The un- 
certainty on the tracking efficiency is given by the statistical precision and the knowledge 
of the purity of the sample oi J/ip ^ /i/i candidates. The purity is estimated by fitting 
the di-muon invariant mass distribution of the J/ijj /j/j candidates with a Crystal Ball 
function [26] to describe the signal shape and a linear background. An alternative estimate 
is obtained by fitting only the linear background on either side of the di-muon resonance. 
The difference in the efficiency evaluated using the two purity methods is taken as the 
systematic uncertainty. 

All particles pass through approximately 20% of a hadronic interaction length of ma- 
terial prior to the final tracking station. Early showering of hadrons reduces the hadron 
tracking efficiency compared to the muon tracking efficiency. An additional correction 
factor to the muon tracking efficiency of (84.3 ± 1.5)% for hadrons is applied which has 
been estimated using the full detector simulation, where the uncertainty on this correction 
corresponds to an uncertainty of 10% in the material budget [27]. 

The electron tracking efficiency uses a tag-and-probe method on Z — > ee data events. 
The tag must satisfy all the electron reconstruction and identification requirements and 
the probe is selected as the highest energy ECAL cluster in the event not associated with 
the tag. The purity of the sample is found, from simulation, to depend on the px of 
the tag. The dependence of the purity is fitted with signal and background templates 
obtained from same-sign and opposite-sign events from data. No momentum information 
is available for the probe, so the tag-and-probe technique only provides an overall tracking 
efficiency for the electrons, which is measured to be (83±3)%. The momentum dependence 
is taken from Z — >■ ee and Z — >■ rr simulation. The electron tracking efficiency uncertainty 
is taken from the fit uncertainty added in quadrature to the statistical uncertainty. 

The identification efficiency, Eid, is measured for muons with the tag-and-probe method 
for the Z —i- fifi data, using a reconstructed track as the probe lepton and evaluated as 
a function of the probe momentum. For low momenta the efficiency is evaluated using 
a, J/ip ^ /i/i sample as a function of both probe pseudorapidity and momentum. The 
muon identification efficiency is found to vary between 93% and 99% in pseudorapidity 
and momentum. The muon Sid uncertainty is evaluated with the same method used for 
the muon £trk uncertainty. 

The electron identification efficiency is measured as a function of probe momentum 
using the tag-and-probe method on Z — >■ ee data and simulation events. The isolation 
requirement introduces a bias of 1% — 4% in data and reconstructed simulation and 
so simulation without the isolation criteria is used instead. The electron identification 
efficiency is found to vary between 85% and 96%, with an uncertainty in each bin estimated 
as the difference in the biased efficiencies from data and simulation. 

The hadron identification efficiency is determined using events triggered on a single 
VELO track. The highest track in each minimum bias event is assumed to be a 
hadron, as verified by simulation. The hadron identification efficiency is taken as the 
fraction of tracks fulfilling the hadron identification requirements. Although the minimum 
bias topology differs significantly from the Z — )■ rr topology, an efficiency dependence 
is observed only in pseudorapidity and so the efficiency is evaluated as a function of 
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pseudorapidity and found to vary between 92% and 95%. The uncertainty for each bin of 
pseudorapidity is estimated as the statistical uncertainty of the bin. A summary of the 
systematic uncertainties is given in Table 2. 



6 Results 

The cross-sections for each analysis stream are determined using Eq. 1, the values given 
in Table 1, and the systematic uncertainties presented in Table 2. The results are 

Upp^z-^rr (r^T^) = 77.4 ± 10.4 ± 8.6 ± 2.7 pb 

Gpp^Z^rr (r^Te) = 75.2 ± 5.4 ± 4.1 ± 2.6 pb 

Gj^^Z^rr {TeT^) = 64.2 ± 8.2 ± 4.9 ± 2.2 pb 

app^z-,rr {r^Th) = 68.3 ± 7.0 ± 2.6 ± 2.4 pb 

a^^Z^rr {TeTh) = 77.9 ± 12.2 ± 6.1 ± 2.7 pb 

where the first uncertainty is statistical, the second uncertainty is systematic, and the 
third is due to the uncertainty on the integrated luminosity. 

A global fit is performed using a best linear unbiased estimator [28] including correla- 
tions between the final states, and a combined result of 

(Tpp^z^TT = 71.4 ± 3.5 ± 2.8 ± 2.5 pb 

is obtained, with a P^r degree of freedom of 0.43. The statistical uncertainties are 
assumed to be uncorrelated as each analysis stream contains mutually exclusive datasets. 
The luminosity and any shared selection or reconstruction efficiencies are assumed to be 
fully correlated. 

A graphical summary of the individual final state measurements, the combined mea- 
surement, the Z — ^ /iyU measurement of Ref. [5], and a theory prediction is shown in Fig. 2. 
The theory calculation uses Dynnlo [29] with the MSTW08 next-to-next-leading-order 
(NNLO) PDF set [23], and is found to be 74.3l^;? pb. 

The ratio of the combined cross-section to the LHCb Z ^ jiji cross-section measure- 
ment [5] is found to be 

where the uncertainty is the combined statistical, systematic, and luminosity uncertainties 
from both measured cross-sections, which are assumed to be uncorrelated. 



7 Conclusions 

Measurements of inclusive Z — )■ rr production in pp collisions at ^/s = 7 TeV in final states 
containing two muons, a muon and an electron, a muon and a hadron, and an electron and 
a hadron have been performed using a dataset corresponding to an integrated luminosity 
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Table 2: Systematic uncertainties expressed as a percentage of the cross-section for each 
Z ^ TT analysis stream. Contributions from acceptance A, branching fractions i3, num- 
ber of background events A^bkg, reconstruction efficiencies Srec; and selection efficiencies 
£sei are listed. The superscripts on £trk^*^ and £id^*^ indicate the first or second r lepton 
decay product candidate. The percentage uncertainties on the cross-section for A^bkg are 
quoted for each individual background, as well as the total background. The efficiency 
uncertainties are split in a similar fashion. 



Stream 






[%] 












' e ' ^ 






A 


1 48 


1.61 


1.32 


1 1 n 


1 1 1 
1.11 


B 


0.46 


0.32 


0.32 


0.32 


0.33 




-^QCD 


4.33 


80 


3 08 

O . WO 


0.40 


0.92 




-^EWK 


4.22 


1.54 


1.52 


0.40 


0.72 


iVbkg 


Ntt 


0.02 


0.08 


0.12 


0.00 


0.58 




Nww 


0.02 


0.14 


0.13 


0.09 


0.08 




Nz 


8.00 






0.22 


0.23 


Total A^bkg 


10.03 


1.75 


3.44 


0.61 


1.32 




£^GEC 


0.10 


0.10 


0.10 


0.10 


0.10 




etrg 


0.88 


0.71 


2.29 


0.72 


4.30 




etrk(l) 


0.71 


0.74 


3.67 


0.79 


3.67 


£^rec 




0.34 


3.67 


0.61 


1.76 


1.68 






0.38 


0.28 


1.72 


0.29 


1.73 






0.78 


0.18 


0.56 


0.03 


0.09 


Total 


^rec 


1.47 


4.21 


4.73 


2.08 


6.15 




^kin 




1.04 


2.89 




1.91 






1.79 


1.91 


3.19 


1.65 


2.75 






1.08 


1.03 


1.86 


0.60 


0.97 




eips 


2.70 






1.92 


2.85 






2.03 










Total 


^sel 


3.97 


2.41 


4.69 


2.60 


4.50 


Total systematic 


11.13 


5.41 


7.56 


3.88 


7.88 
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Figure 2: Measured cross-sections for the Z decaying to the final states t^t^, T^T"e, rgX^, 
r^T/i, and TeT/i (open points) compared with theory (yellow band) and the combined Z — )• 
TT and LHCb Z — >• /i/i measurements (closed points) where pt and ?7 are the transverse 
momentum and pseudorapidity of the leptons, and M is the di-lepton invariant mass. 
The inner error bars represent statistical uncertainty while the outer error bars represent 
combined statistical, systematic, and luminosity uncertainties. The central theory value 
is given by the light yellow line while the associated uncertainty by the orange band. 
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of 1028 ± 36 pb . The cross-sections for the individual states have been measured in 
the forward region of 2.0 < t]^ < 4.5 with > 20 GeV/c and 60 < M^r < 120 GeV/c^, 
and a combined result calculated. The results have been compared to Standard Model 
NNLO theory predictions and with the LHCb Z — )> yUyU cross-section measurement. The 
individual measurements, the combined result, the Z — > fi^ cross-section, and the theory 
prediction are all in good agreement. The ratio of the Z ^ iifi cross-section to the Z ^ tt 
cross-section is found to be consistent with lepton universality. 
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